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We show that the origin of electronic transitions of molecular many-body systems can be revealed 
by a quantified natural transition orbitals (QNTO) analysis and the electronic excitations of the 
total system can be mapped onto a standard orbitals set of a reference system. We further illustrate 
QNTO on molecular systems by studying the origin of electronic transitions of DNA moiety, thymine 
and thymidine. This QNTO analysis also allows us to assess the performance of various functionals 
used in time-dependent density functional response theory. 
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Quantum mechanical description of electronic excita- 
tions of molecular many-particle systems has been an im- 
portant topic to date. Time-dependent linear response 
DFT (TDDFT) 0, Q has been widely used to calculate 
the electronic excitation energies and oscillator strengths 
of large molecular systems. In contrast, high-level quan- 
tum chemistry approaches such as equation of motion 
coupled cluster singles and doubles (EOM-CCSD) Q can 
deal with only small systems. The TDDFT results can 
be, however, unreliable due to the quality of exchange- 
correlation functionals used. While assessment of differ- 
ent functionals often focuses on the excitation energies 
and oscillator strengths 0, [f| , verification of the phys- 
ical origin of electronic transitions should be desirable. 
Casida has proposed that TDDFT wavefunction can be 
interpreted in terms of linear combination of singly ex- 
cited configurations (LCSEC) Q. Casida's scheme sug- 
gests that a systematic analysis can be developed so that 
LCSEC of TDDFT and many-body theories like EOM- 
CCSD can be linked for each excitation. The LCSEC 
contains information on transition origin that directly 
connects to the electronic structure change which governs 
the physics and chemical actions, like nuclear dynam- 
ics and chemical reactions, of an electronically excited 
molecule. 

For assigning the transition origin of an excitation from 
LCSEC, the most important property is its universality. 
In other words, each excitation can only have one unique 
representation of transition origin. In this regard, we 
find that natural transition orbital (NTO) analysis [6] is 
a good starting point because it relies on the minimum 
representation of the transition operator using a NTO 
pair consisting of an electron- and a hole-orbital for an 
excitation. However, we also find that in many molecular 
systems, an excitation cannot be described very well by 
single NTO pair. One of the typical and simplest exam- 
ples is a well separated dimer. Therefore, we develop a 
quantified NTO (QNTO) analysis that allows taking into 



account multi-NTO pairs properly as needed. Moreover, 
NTO analysis uses a pictorial description of a NTO pair 
that can result in ambiguity especially for large molecules 
with no symmetry. To establish a clear interpretation of 
transition origin of NTO pairs, we introduce a standard- 
orbitals set onto which the hole- and electron-orbital of 
a NTO pair can be mapped. This scheme can be applied 
for a wide range of molecular systems, such as DNA and 
proteins. With QNTO analysis, we demonstrate that 
separate dimer can be properly described by two NTO 
pairs. In addition, the variation of the transition origin 
of an excitation due to the environment change of the 
system can also be determined based on the coefficient 
change of standard-orbitals. Furthermore, we show that 
TDDFT together with recently developed long-range cor- 
rected (LC) hybrid functionals 0,0, &E1 (LC-TDDFT) 
can predict electronic transition origins of similar quality 
as EOM-CCSD. 

For single-reference linear response methods for elec- 
tronic excitations, e.g., configuration interaction singles 
(CIS), time-dependent Hartree-Fock (TD-HF), EOM- 
CC, and TDDFT, the NTO analysis @ is known to give a 
concise transition picture out of the original many SECs 
of an (often low-lying) excitation. This is due to the fact 
that mathematically some SECs in a LCSEC can be com- 
bined to form the minimum representation of transition 
operator. For example, let us consider an electronic exci- 
tation whose LCSEC is 0.3*^ + 0.5*? + 0.8*^ where, 
for example, *™ = a^Si* and * is the ground state 
wavefunction of a theoretical method which may be a 
single determinant (e.g. in DFT) or a combination of 
multi-determinants (e.g. in EOM-CCSD). If one follows 
a NTO-generating procedure - a unitary transformation 
to the original occupied and virtual orbitals - the best re- 
duced LCSEC can be achieved with an additional phase 
confirmation of SECs. For the above example, it reduces 
to 0.3*? + 0.5*| + 0.8*5 = V0.9172*^ + V0.0628*f 
where ai = Ni(cx&i + c^ay), aj = Nj(—C2(ii + C\a,j), 
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a\ = N A{d\a\+ d 2 a\) and &jg = N B {— d 2 al + dial) with 
ci = 0.5696, c 2 = 0.8219, rfi = .1784, d 2 = 0.9840, 
and, for instance, Nj = 1/ y^cf + cf . In this new descrip- 
tion, the first SEC tyf, corresponding to the first dom- 
inant NTO pair (NTOl), now dominates and therefore 
the transition picture of the excitation is clear. 

However, as above mentioned, there may be a case in 
which the NTOl cannot be representative to an excita- 
tion, e.g. in a separate dimer. For this reason, in the 
first step, we generalize the traditional NTO analysis by 
taking into account the 2nd dominant NTO pair (NT02) 
if the NTOl has less than 70% domination of total LC- 
SEC. We also establish a procedure to determine the rel- 
ative phase between NTOl and NT02. Adding NT03 or 
higher term can readily be done. For the physical prop- 
erties associated with one-particle operators like elec- 
tronic density, the phase difference leads to no difference, 
(tt? ± ^\p(r)m ± *§} = (*?|#r)|*?> + (*5|p(r)|*5), 
i =/= j, while it is not for two-particle operators such as 
the interelectron coulomb repulsion. 

As the 2nd step of generalization, we project the hole- 
(NTOl-H) and electron-orbital (NTOl-E) of NTOl onto 
a chosen standard orbitals set to see their quantitative 
involvement. The same procedure will be performed for 
NT02. In doing so three advantages emerge: 

(1) The precise contribution of a particular orbital- 
type of interest in an excitation can be computed. 

We consider 0.3*f+0.5*"+0.8*J as an example whose 
NTOl is fyf, if we simply use {ipi 7 ipj 7 ip a ,ipi,} as the 
standard-orbitals, the proportions of ifi and (fj in the 
NTOl-H are 0.1784 and 0.9840, respectively Similarly, 
the proportions of ip a and tp b in the NTOl-E are 0.5696 
and 0.8219, respectively Having these results, say, if ipi 
and ipj originate from the holc-orbitals of 1 nir* and 1 tttt* 
transitions, respectively, then the precise proportion of 
the two orbitals in the excitation of interest can be ob- 
tained. Since NTO-based LCSEC is unique, the orbital 
projection would not be misleading. 

(2) Improving the pictorial NTO interpretation. 
Although NTO analysis can largely reduce many SECs 

in the original LCSEC to the dominant one or two, it 
is only rendered as an assistant to see the qualitative 
transition picture by plotting NTOl-H and NTOl-E. On 
the other hand, projecting NT01(2)-H(E) to standard- 
orbitals set would provide a rigorous explanation of tran- 
sition origin. 

(3) Unification of orbital-standard for comparison. 

If the exact ground state wavefunction of an electronic 
system can be approximated by either \^q S ) of DFT or 
of CCSDj the wavefunction up to LCSEC for 
an electronic excitation can be obtained within the linear 
response theory. However, the LCSECs given by TDDFT 
and EOM-CCSD arc based on different references, thus 
making a direct comparison difficult. Furthermore, am- 
biguity can also arise if one compares LCSECs for excita- 
tions in two systems with different environments. Thcrc- 



TABLE I. Recovery of exciton theory prediction (TD- 
a;B97X) for (a) 1 ^^* and (b^mr* excitations. 'P' denotes 
the phase of NT02. In (a) as the C6-C5' distance, defined 
in JTfJj increases, the two local transitions (NTOl, NT02) 
tend to have a similar domination. The phase flipping be- 
tween 9.24A and 10.24A may be due to a conical intersection. 
In contrast, in (b) two local transitions remain uncoupled re- 
gardless of the C6-C5' distance. 



(a) 


NT01:"P1„ - SV 




NT01:"P1 U - SV 




C6-C5' 


NT02:"P1 U - Sl„" 


P 


NT02:"P1„ - SV 


P 


(A) 


A f 1st 2nd 
(nm) (%) (%) 




A f 1st 2nd 
(nm) (%) (%) 




4.24 


230.4 0.104 65 32 




228.4 0.227 65 31 


+ 


5.24 


229.6 0.114 69 28 




228.4 0.242 69 28 


+ 


6.24 


229.3 0.118 68 28 




228.5 0.252 68 28 


+ 


7.24 


229.2 0.119 65 32 




228.7 0.261 65 32 


+ 


8.24 


229.1 0.118 57 39 




228.7 0.268 57 39 


+ 


9.24 


229.0 0.118 50 46 




228.8 0.272 50 46 


+ 


10.24 


228.8 0.270 64 33 


+ 


229.0 0.123 64 33 





(b) 
C6-C5' 

(A) 


NT01:"N1„ - Sl u " 

A f 1st 
(nm) (%) 


NT01:"N1„ - SV 
A f 1st 
(nm) (%) 


4.24 


243.9 0.0001 99 


242.14 0.0002 99 


7.24 


242.8 0.0002 100 


242.10 0.0002 100 


10.24 


242.5 0.0002 100 


242.11 0.0002 100 



fore, importantly, projection of NT01(2)-H(E) of LCSEC 
onto a standard-orbitals would make it possible to com- 
pare two LCSECs based on different references. 

We demonstrate that QNTO analysis correctly recov- 
ers the molecular exciton theory [hi ] prediction in which 
a molecular homo dimer system composed of well sepa- 
rated geometrically identical two monomers (ideal dimer) 
has the formation of two excited states, -^(ft'fiv ± 

ip u ipl), where ip^ and <p u is the excited and unexcited 
wavefunctions of monomer u and ip\ and <p v those of 
monomer v. This situation should be reflected in QNTO 
analysis by having NTOl and NT02 as each monomer's 
local transition with close domination. QNTO helps ver- 
ify if a real dimer system can be appropriately described 
by the exciton model or not. Table I shows this trend of 
two 1 irir* excitations and no similar trend for two 1 nw* 
ones of a real homo dimer system composed of two ge- 
ometrically slightly different monomers in various sepa- 
rations. In addition, QNTO can be applied to molecular 
hetero dimer as well. 

After generalizing NTO to QNTO analysis, obtaining 
a reliable LCSEC for an excited state is the next im- 
portant concern. It depends on the use of proper ap- 
proximate exchange-correlation functionals in TDDFT 
calculations mentioned earlier. In this regard, recently 
developed LC hybrid functionals, are proved to be good 
candidates to describe long-range transitions, and some- 
times even outperform the traditional (non-LC) hybrid 
functionals for short-range transition as well. In a typi- 
cal LC hybrid functional, 100% Hartree-Fock exchange 
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FIG. 1. Molecular structrues of (a) single thymine (Thy) and 
(b) thymidine (dT) extracted from ideal B-DNA. The dT has 
an additional backbone structure of sugar connecting to Thy. 



(a) 



PI 



Nl P2 



N2 



SI S2 



FIG. 2. Standard-orbitals generated in the Thy B3LYP-DFT 
calculation. PI, Nl, P2 and N2 correspond to the HOMO, 
HOMO-1, HOMO-2, and HOMO-3, respectively, and SI and 
S2 are LUMO and LUMO+1. PI and P2 (SI and S2) are 
the major components of the h- (e-)orbital of 1 7T7r*-character 
excitations, while Nl and N2 are those of h-orbital of rm*- 
character excitations. 



is employed for a long-range part of the interclcctron 
repulsive operator L(r)/r, while a GGA exchange func- 
tional is employed for the complementary short-range 
part. Recently, Chai and Head-Gordon found that LC 
hybrid functionals can outperform global hybrid func- 
tionals (e.g. B3LYP) in several important applications, 
such as thermochemistry, thermal kinetics, noncovalent 
systems, dissociation of symmetric radical cations and 
long-range charge-transfer excitations between two well- 
separated molecules Recent TDDFT performance 
tests of various LC- and non-LC hybrid functionals for 
predicting excitation properties of several molecular sys- 
tems have been reported. Jacqucmin et al. assessed the 
accuracy of wB97 family for predicting excitation ener- 
gies [i| , whereas Caricato et al. focused on the simulated 
spectra of three molecular groups - Alkenes, Carbonyls 
and Azabenzenes - in reference to EOM-CCSD results 
[Bj]. Since here we emphasize NT01(2) transition origin 
of electronic excitations, several excitations of two bench- 
mark systems - thymine (Thy) and thymidine (dT) - cal- 
culated by various theoretical methods are investigated 
as examples. 

The prepared Thy and dT (Fig. 1) are extracted from 
ideal B-DNA structure. EOM-CCSD and TDDFT sin- 
glet excitations calculations are carried out using Gaus- 
sian09 package |15( with 6-31G(d) basis set. For simplic- 
ity, solvent models are not used. EOM-CCSD and TD- 
wB97X Q are demonstrated here whereas several other 
theoretical methods are tested in . 

The results of several chosen excitations of Thy and dT 
are shown in Fig. 3, Table II and [l(|. The chosen exci- 
tation, regarded as SO-hosted (standard-orbitals hosted), 
has its NTOl-H and NTOl-E overall more than 30% den- 
sity contributions from the monomer orbitals (standard- 
orbitals) shown in Fig. 2(a) and (b), respectively. To 



TABLE II. Electronic transition origin (shown in the NTOl 
column) of the first 7 SO-hosted singlet excitations of ideal 
Thy and dT given by EOM-CCSD and TD-ojB97X calcula- 
tion. N denotes the excitation order. % records the ratio of 
NTOl to the excitation. Type denotes the excitation clas- 
sification based on EOM-CCSD Thy NTOl expressions; for 
EOM-CCSD the total domination of LCSEC (%) is also listed 
in the same column. 
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NTOl 


% 


Type 


N|NT01 


\%\ 


Type 


EOM- 
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Nl - SI 


85 


A(86) 


1 Nl - SI 


79 


A(80) 


CCSD 
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PI - SI 


83 


B(83) 


2 PI - SI 


83 


B(83) 
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Nl N2 - 


86 


C(88) 


3 Nl N2 - 


72 


C(75) 






(SI) S2 






(SI) S2 
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P2 - SI 


83 


D(83) 


4 P2 - SI 


69 


D(83) 
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PI - S2 


86 


E(86) 


5 PI - S2 


75 


E(84) 
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N2 - SI 


83 


F(85) 


6 N2 (B) - SI 


78 


(F)(81) 
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Nl (N2) 


72 


G(84) 


7 Nl (N2) - 


71 (G)(75) 






- S2 






(SI) S2 






TD- 


1 


Nl - SI 


100 


A 


1 Nl - SI (S2) 99 


(A) 


wB97X 


2 


PI - SI 


97 


B 


2 PI - SI 


97 


B 




3 


Nl N2 - 


99 


C 


3 Nl N2 - 


98 


C 






(SI) S2 






(SI) S2 
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P2 - SI 


99 
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4 PI - S2 


96 


E 
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PI - S2 


98 
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5 P2 - SI 


97 


D 
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N2 - SI 


97 
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6 N2 B (O) 


96 


(F) 












- SI 
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Nl N2 - 


91 


(G) 


7 Nl N2 - 


83 


(G) 






(SI) S2 






(SI) S2 







confirm this, their projection to standard-orbitals is first 
performed and the density contribution of a standard- 
orbital is given by the square of its coefficient; a transition 
expression is also determined from the coefficients (de- 
tailed in [lH). When dealing with dT, wc also consider 
the involvement of all orbitals originated from the back- 
bone structure (denoted as 'B') to a NT01(2)-H. The 
integration of standard-orbitals (plus 'B' for dT) gen- 
erally accounts for more than 90% of a NT01(2)-H(E). 
Otherwise, we attribute the involvement of all other Thy- 
orbitals as 'O'. Moreover, a precise statistical assessment 
of NTOl of an excitation is carried out as well, where 
the root mean square deviation of the projection coef- 
ficients (a) of NTOl-H and NTOl-E of each electronic 
transition is computed, om is used when we map the 
NTOl obtained from TDDFT with various functionals 
onto NTOl of EOM-CCSD, and a E is for comparison 
between EOM-CCSD results for Thy and dT. 

From Table II, it can firstly be seen in EOM-CCSD 
results that the low-lying excitations indeed can be well 
described by LCSEC, as its total domination to whole 
excitation is generally higher than 0.80 (0.75 to be the 
lowest); secondly their LCSEC can generally be well ex- 
plained by NTOl, as its domination is often close to the 
LCSEC one (NT02 may arise for higher-lying excitations 
or the dimer systems mentioned earlier). TD-wB97X 
(and other LC-TDDFT) results agree quite well with 
those of EOM-CCSD, even for the small system Thy, not 
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(b) 



(c) 



^0.2 ■ 

w 
O 



- EOM-CCSD 
TD-(oB97X 
TD-B3LYP 
TD-LSDA 
TD-HF 



dT 




TD-iuB97X 
TD-B3LYP 
TD-HF 



dT 



A B C D E F G 

Thy 



A B C D E F G 

* .Thy 



180 200 220 240 

Wavelength (nm) 



A B C D E F G 

Excitation Type 



A B C D E F G 
Excitation Type 



FIG. 3. (a) Absorption wavelengths calculated using TD- 
LDA (pure functional), TD-HF (exact exchange), TD-B3LYP 
(global hybrid functional) and TD-o;B97X (LC hybrid func- 
tional) vs. absorption wavelengths from EOM-CCSD. Data 
of each method for the same excitation lies on the same verti- 
cal line, (b) Oscillator Strength: line denotes the results from 
EOM-CCSD calculation, (c) a M of the Type-A~G calculated 
by the same methods except TD-LDA; erg is also plotted along 
with aM in the Thy result. Type-G is not predicted by TD- 
HF calculation for Thy whereas the excitations of dT given 
by TD-LDA have quite different transition origin and cannot 
be well-classified as Type-A~Type-G [l(|. 



only in terms of the absorption energy, oscillator strength 
but also in terms of the predicted NTOl transition ori- 
gin (Fig. 3), and thus have notably outperformed other 
tested theoretical methods [HI]. Incidentally, one should 
note that EOM-CCSD provides generally smaller NTOl 
domination than TDDFT does even for the lowest exci- 
tation of Thy. This is attributed to the limitation of adi- 
abatic approximation of exchange-correlation kernel used 



in this work [17|. In other words, TDDFT renders LC- 
SEC to be 100% of an excitation. Furthermore, LCSEC 
from LC-TDDFT agrees well with that of EOM-CCSD 
in that the relative percentage of different SECs is well 
reproduced. Therefore, the a M of LC-TDDFT for the 
low energy absorptions in Table II is very small. 

The key factor causing functionals to have the re- 
sults of different quality is their HF-exchange prescrip- 
tion, as shown in [l6| . Global hybrid functionals, like 
B3LYP and PBE1PBE, bear increased self-interaction er- 
ror (SIE) for charge-transfer transition component with 
systematic underestimation of transition energy. Instead, 
LC hybrid functionals, implemented with a distance- 
dependent HF-exchange, can efficiently reduce SIE. In 
ll|, it has been argued that LDA and GGAs (semilocal 
density functionals) are accurate only for small interelec- 
tronic separations (r). Even though how precisely the 
amount of HF-exchange should change with r remains 
unknown (l3j , the importance of HF-exchange generally 
grows with an increasing r. Compared with the results of 



other theoretical methods demonstrated in [16[ we find 
that TD-wB97X, having the smallest average <tm, best 
reproduces the NTOl transition origin that EOM-CCSD 
predicts for the several low-lying excitations of Thy and 
dT without sacrificing other excitation properties. Hence 
it is a very good choice to use when studying the elec- 
tronic excitations of DNA systems. For other molecular 
systems more tests may still be needed; for open-shell 
systems the use of present model also requires extra con- 
sideration where the ground state wavefunction is often 
no-longer dominated by a single Slater determinant. 

For electronic excitations, QNTO analysis provides a 
definite assessment of their transition origin which has 
direct connection to the electronic structure upon which 
later molecular dynamics (nuclear motions, non-radiative 
transitions, photo-induced chemical reactions, and etc.) 
depend. When one deals with several electronic excita- 
tions of a molecule with low or even no symmetry, it also 
provides the fingerprint for them. A measure of transi- 
tion origin variation of particular electronic excitations 
(e.g. SO-hosted ones in this work) due to different im- 
plemented theoretical methods (cm) or different environ- 
ments (<te) can be performed as well. Since the compu- 
tational cost is moderate, LC-TDDFT cooperated with 
QNTO analysis is expected to be particularly useful in 
studying larger systems. 

We thank National Science Council and NCTS of Tai- 
wan for supports. 
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Supplemental Material: The Quantified NTO Analysis for the Electronic Excitations 

of Molecular Many-Body Systems 




(a) 



FIG. SI. The Thy dimer used in Table I for the demonstration 
of molecular exciton theory recovery by quantified natural 
transition orbitals (QNTO) analysis. In the process the dis- 
tance of C6-C5' (dash line) is modified with the two Thy base 
planes bearing unchanged polar angle and azimuthal angle in 
spherical coordinate system of z-axis along C6-C5'. In other 
words, the transition-dipole-transition-dipole-product of an 
excitation, a quantity on which the exciton theory depends, is 
fixed and only the distance between the two transition-dipoles 
is differentiated. 
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- CAM-B3LYP 
■PBE1PBE 

■ B3LYP 




r (Bohr) 



FIG. S2. L(r) vs. r for various functionals. L(r), defined as 
(HF operator * r), is a fraction of Hartree-Fock exchange at 
interelectronic separation r. 
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FIG. S3. The results of more theoretical methods 
not mentioned in Fig. 3. Notations are detailed in 
the main text. (a)Absorption wavelengths calculated us- 
ing TD-CAM-B3LYP(») (U, TD-cjB97X-D(.) [H, TD- 
cjB97(») [13, TD-LC-^PBE(») [13 (LC hybrid function- 
als), TD-PBEIPBE(A) [H (a global hybrid functional) and 
CIS(*) Hi] (TD-HF with Tamm-Dancoff approximation) vs. 
absorption wavelengths from EOM-CCSD(/) [gj. Data of 
each method for the same excitation lies on the same vertical 
line, (b) Oscillator Strength: line denotes the results from 
EOM-CCSD calculation, (c) a M of the Type-A~Type-G ex- 
citations. 
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TABLE SI. Electronic transition origin of the first 8 SO- 
hosted singlet excitations of ideal Thy and dT given by calcu- 
lation of (a) TDDFT with four different LC hybrid function- 
als: CAM-B3LYP, ljB97, wB97X-D and LC-wPBE and (b) 
other theoretical methods: CIS, TD-HF, TD-LDA [H, |S^ 
and TD-PBE1PBE not demonstrated in Table II. Inciden- 
tally, in CIS result of dT only 7 SO-hosted excitations out of 
the first 50 states (to 99.14nm for the 50th one) are recog- 
nized and listed. N and P, respectively, denote the excitation 
order and the phase of NT02 for the absorptions with NTOl 
domination to the linear combination of singly excited con- 
figuration (LCSEC) lower than 70%. The next column shows 
the transition origin of NT01(2), followed by the NTOI(2) 
domination to the excitation. The last column denotes the 



excitation classification based on EOM-CCSD Thy NTOl ex- 
pressions. The excitation denoted as (F') has similar Type-F 
NTOl expression but larger om- 
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Ideal dT 




N; 


NTOl 


% 


Type 


N; 


NT01(2) 


% 


Type 




P 








P 








TD- 


1 


Nl - SI 


100 


A 


1 


Nl - SI 


99 


A 


CAM- 


2 


PI - SI 


97 


B 


2 


PI - SI 


97 


B 


B3LYP 


3 


Nl N2 - (SI) S2 


98 


C 


3 


Nl N2 - (SI) S2 


97 


C 




4 


P2 - SI 


95 


D 


4 


PI - S2 


09 


E 




5 


PI - S2 


94 


E 


- 


P2 - SI 


29 






6 


N2 - SI (S2) 


95 


(F) 


5 


P2 - SI 


71 


D 




7 


Nl N2 - S2 


94 


(G) 


6 


N2 B - SI (S2) 


92 


(F) 




9 


P2 - S2 


93 


- 


7 


N2 B - SI 


91 


(F') 












8 


Nl N2 - S2 


86 


(G) 


TD- 


1 


Nl - SI 


100 


A 


1 


Nl - SI 


99 


A 


uB97X 


2 


PI - SI 


96 


B 


2 


PI - SI 


97 


B 


-D 


3 


Nl N2 - (SI) S2 


97 


C 


3 


Nl N2 - (SI) S2 


96 


C 




4 


P2 - SI 


94 


D 


4 


PI (P2) - (SI) S2 


56 


(E) 




5 


PI - S2 


93 


E 


- 


(PI) P2 - SI (S2) 


43 






6 


N2 - SI (S2) 


94 


(F) 


5 


(PI) P2 - SI (S2) 


59 


(D) 




7 


Nl N2 - S2 


94 


(G) 


+ 


PI (P2) - (SI) S2 


39 






9 


P2 - S2 


94 




6 


N2 B - SI (S2) 


91 


(F) 












7 


N2 B - SI 


89 


(F') 












8 


Nl N2 - S2 


X:i 


(G) 


TD- 


1 


Nl - SI (S2) 


99 


(A) 


1 


Nl - SI (S2) 


99 


(A) 


uB97 


2 


PI - SI 


97 


B 


2 


PI - SI 


97 


B 




3 


Nl N2 - S2 


99 


(C) 


3 


Nl N2 - S2 


98 


(C) 




4 


P2 - SI 


94 


D 


4 


PI - S2 


90 


E 




5 


PI - S2 


93 


E 


5 


P2 - SI 


91 


D 




6 


N2 (O) - SI 


94 


(F) 


6 


N2 (B) (O) - SI 


90 


(F) 




7 


Nl (N2) - (SI) S2 


85 


(G) 


7 


Nl (N2) - (SI) S2 


81 


(G) 




9 


P2 - S2 


94 




8 


(N2) B (O) - SI 


80 


(F') 


TD-LC- 


1 


Nl - SI (S2) 


1110 


(A) 


1 


Nl - SI (S2) 


99 


(A) 


uPBE 


2 


PI - SI 


97 


B 


2 


PI - SI 


97 


B 




3 


(Nl) N2 - S2 


99 


(C) 


3 


(Nl) N2 - S2 


99 


(C) 




4 


P2 - SI 


94 


D 


4 


PI - S2 


92 


E 




5 


PI - S2 


92 


E 


5 


P2 - SI 


93 


D 




6 


N2 (O) - SI 


95 


(F) 


6 


N2 (B) (O) - SI 


97 


(F) 




7 


Nl (N2) - (SI) S2 


85 


(G) 


7 


Nl (N2) - (SI) S2 


83 


(G) 




8 


P2 - S2 


94 




8 


(N2) B (O) - SI 


78 


(F') 



Note for Table SI: 



The coefficient threshold to output singlet excited con- 
figurations (SECs) for QNTO analysis in TDDFT calcu- 
lations is set to 0.01; for EOM-CCSD it is set to 0.1. 
Even though the threshold used is different, we assume 
the SECs with coefficient smaller than 0.1 have a minor 
effect on the composition of LCSEC so that EOM-CCSD 
results can still be used as a norm for TDDFT results 
for comparison. In Table SI NT01(2) is the acronym of 
the (2nd) dominant natural transition orbital pair. PI, 
Nl and etc. are the standard-orbitals defined in Fig. 
2 used to expand the hole- (NT01(2)-H) and electron- 
orbital (NT01(2)-E) of a NT01(2). To determine the 
transition expression of a NT01(2), PI, for example, is 
involved, partial-involved and not involved in a NTOl- 
H if the \(P1\NT01 - iJ)| 2 , density contribution of PI 
to the NTOl-H, is [1,0.3], (0.3,0.1] and (0.1,0] respec- 
tively; (P\\NTO\ — H) denotes the projection coefficient 
of the NTOl-H onto PI. Based on these coefficients 
the SO-hosted excitations, assigned as having NTOl- 
H and NTOl-E overall more than 30% density contri- 
butions from the standard-orbitals shown in Fig. 2(a) 
and Fig. 2(b), respectively, can also be readily identi- 
fied from several excitations calculated by a theoretical 
method. Inside " ()" in a transition expression indicates 
the standard-orbital that has partial-involvement; the 
underline below a standard-orbital means its coefficient 
is negative. Moreover, when dealing with dT we also in- 
troduce 'B' referring to all orbital- fraction of NTOl-H 
from the DNA sugar backbone attaching to Thy. If the 
integrated density contribution of standard-orbitals (plus 
'B' for dT) is less than 90% to a NT01(2)-H(E), the re- 
maining part that is from other Thy-orbitals is denoted 
as 'O' in the transition expression (shown in "()" if its 
contribution is (0.3,0.1]). 

It also requires extra cautions if the comparison of 
two transition expressions with some phase difference of 
standard-orbitals is being carried out. For instance, pay- 
ing attention on the TD-PBE1PBE result, one finds that 
the NTOl represented by "Nl N2 - (SI) S2" for the 5th 
excitation of Thy is not fully equal to the NTOl, "Nl 
(N2) - (SI) S2" , for the 5th excitation of dT because 
of the SI phase difference; even single-particle operator, 
e.g. /5(r), can recognize the difference. It can be proved 
that (*? + *f|p(r)|*? + *, 6 ) + m - *J|p(r)|*? - ttj) 
because (^f |a^a.t,|\E , |') connects the cross terms. On the 
other hand, the relative phase between NTOl and NT02 
can be absorbed into NTOl or NT02. For example, 
the complete expression for the 7th excitation of dT in 
the TD-PBE1PBE result reads VolM"(Pl) N2 (B) - SI 
(S2)"- V / 0T35"P1 (Nl) B - (SI) S2". Absorbing the rela- 
tive phase between NTOl and NT02 into NT02-H, for 
example, it can be rewritten as V0.64" (PI) N2 (B) - SI 
(S2)"+V035"P1 (Nl) B - (SI) S2". ('B' phase is not 
taken into account.) 

Results and Discussion 

Among the LC hybrid functionals used in Table SI(a), 
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(b) 


Ideal Thy 


Ideal dT 




N; 


NT01(2) 


% 


Type 


N; 


NT01(2) 


% 


Type 




P 








P 








CIS 


1 


Nl (N2) - SI (S2) (O) 


98 


(A) 


1 


PI - SI 


94 


B 




2 


PI - SI 


94 


B 


2 


Nl (N2) - SI (S2) (O) 


98 


(A) 




3 


(Nl) N2 - S2 (O) 


99 


(C) 


3 


(Nl) N2 - S2 (O) 


97 


(C) 




4 


PI P2 - SI S2 


59 


(E) 


4 


PI - S2 


76 


E 




+ 


PI P2 - SI S2 


36 




5 


P2 - SI (S2) 


81 


(D) 




6 


PI P2 - SI S2 


7:i 


(D) 


7 


P2 - (SI) S2 


87 


_ 




8 


P2 - (SI) S2 


88 


- 


15 


(Nl) N2 O - SI 


47, 


(F) 




12 


Nl N2 (O) - SI 


64 


(F') 


- 


PI (B) - O 


26 






_ 


Nl N2 (O) - S2 (O) 


;i2 














18 


N2 O - SI 


68 


(F) 












+ 


Nl (N2) O - S2 (O) 


23 












TD-HF 


1 


PI - SI 


97 


B 


1 


PI - SI 


97 


B 




2 


Nl (N2) - SI (S2) (O) 


99 


(A) 


2 


Nl (N2) - SI (S2) (O) 


98 


(A) 




3 


(Nl) N2 - S2 (O) 


100 


(C) 


3 


(Nl) N2 - S2 (O) 


98 


(C) 




4 


PI P2 - SI S2 


60 


(D) 


4 


PI - S2 


71 


E 




+ 


PI P2 - SI S2 


36 




5 


P2 - SI S2 


80 


(D) 




5 


PI P2 - (SI) S2 


77 


(E) 


7 


P2 - (SI) S2 


88 


_ 




8 


P2 - (SI) S2 


88 


_ 


15 


(PI) (N2) (B) (O) - SI O 


38 


(F) 




IS 


N2 O - SI 


71 


(F) 


+ 


PI (Nl) (N2) (O) - SI O 


32 






22 


PI O - S2 O 


39 


- 


21 


PI - SI (O) 


82 


- 




+ 


Nl (N2) - O 


32 












TD- 


1 


Nl - SI 


100 


A 


1 


Nl - SI 


100 


A 


B3LYP 


2 


PI - SI 


96 


B 


2 


PI - SI 


97 


B 




3 


N2 - SI 


93 


F 


3 


N2 B - SI 


96 


(F) 




4 


P2 - SI 


92 


D 


4 


P2 - SI 


81 


D 




5 


Nl (N2) - S2 


86 


(C) 


5 


(Nl) (N2) B - SI (S2) 


63 


(C) 




6 


PI - S2 


89 


E 


+ 


Nl B - (SI) S2 


36 






7 


(Nl) N2 - S2 


96 


(G) 


6 


Nl N2 - SI S2 


61 


(CG) 




9 


P2 - S2 


82 


- 


+ 


Nl (N2) B - SI S2 


39 














7 


PI - S2 


75 


E 












9 


(Nl) N2 - S2 


96 


(G) 


TD- 


1 


Nl - SI 


100 


A 


1 


Nl - SI 


100 


A 


PBE1- 


2 


PI - SI 


96 


B 


2 


PI - SI 


1)7 


B 


PBE 


3 


(Nl) N2 - SI (S2) 


!(() 


(F) 


3 


N2 B - SI 


94 


(F) 




4 


P2 - SI 


93 


D 


4 


P2 - SI 


8(1 


D 




5 


Nl N2 - (SI) S2 


80 


(C) 


5 


Nl (N2) - (SI) S2 


(11 


(C) 




6 


PI - S2 


90 


E 


_ 


B - SI (S2) 


37 






7 


(Nl) N2 - S2 


95 


(G) 


6 


PI (Nl) (N2) - (SI) S2 


(,!) 


(E) 




9 


P2 - S2 


90 


- 


+ 


(PI) (N2) P2 (B) - SI (S2) 


30 














7 


(PI) N2 (B) - SI (S2) 


64 


(F) 












- 


PI (Nl) B - (SI) S2 


35 














9 


(Nl) N2 - S2 


1)2 


(G) 


TD- 


1 


Nl - SI 


100 


A 


1 


Nl - SI 


100 


A 


LDA 


2 


PI (P2) - SI 


97 


(B) 


2 


(PI) B - SI 


99 


(B) 




3 


N2 - SI 


100 


F 


3 


PI N2 (B) - SI 


98 


(BF) 




4 


Nl - S2 


98 


(G) 


5 


(PI) N2 B - SI 


!)!) 


(F) 




5 


(PI) P2 - SI (S2) 


92 


(DB) 


6 


Nl - S2 


!)!) 


(G) 




6 


PI - S2 


85 


E 


7 


(PI) P2 - SI (S2) 


85 


(DB) 




7 


N2 - S2 


100 


(C) 


9 


(PI) (N2) B - S2 


86 


(E) 




11 


P2 - S2 


48 




10 


PI (N2) (B) - S2 (O) 


81 


(E) 




+ 


(PI) O - SI 


44 
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LC-cjPBE and wB97 have the same HF-exchange op- 
erator (Fig. S2), while the others have different ones. 
LC-wPBE and o;B97 give quite similar results for both 
Thy and dT, showing that the use of different versions 
of GGA in their functional forms matters little here and 
the fraction of HF-exchange with interelectron distance 
is much more important. We also see that CAM-B3LYP 
and wB97X-D give quite similar results, which has also 
been found in [Sid ] for other molecular systems. 

Fig. S3 and Table SI show that TDDFT with LC hy- 
brid functionals evidently outperform other methods for 
excitation energy, oscillator strength and NTOl tran- 
sition origin prediction. In general, CIS and TD-HF 
severely overestimate the excitation energy and oscilla- 
tor strength, while TD-B3LYP, TD-PBE1PBE and TD- 
LDA underestimate them. The trend that pure density 
functional (LDA) and global hybrid functionals (B3LYP, 
PBE1PBE) tend to give a larger underestimation of os- 
cillator streng th fo r bright absorptions is also observed 
elsewhere of 



Sll 



The values of om for different type 
of excitations yielded by the five methods are in most 
cases larger than those given by TDDFT with LC hybrid 
functionals. Moreover, LDA, having no HF-exchange 
in its functional form, leads to excitations with quali- 
tatively different transition origins that generally cannot 
be well classified by Type-A~Type-G. In fact if we focus 
on Type-B, D and E (bright absorptions), their relative 
absorption energy, oscillator strength, and NTOl transi- 
tion origin can be qualitatively captured by TD-B3LYP 
for both systems. However, for Type-C, F and G dark 
absorptions the description of NTOl transition origin is 
more problematic and for Type-F its excitation energy 
is underestimated more significantly than other absorp- 
tions such that its order goes down; both properties play 
a significant role in later dynamics of an electronically ex- 
cited molecule. TDDFT with LC hybrid functionals, on 
the other hand, can make them (and also excitation en- 
ergy, oscillator strength) agree well with those predicted 
by EOM-CCSD. 

We have observed a close relation between a proper 
description of electronic excitations and the prescrip- 
tion of HF-exchange in different functionals. LC hy- 
brid functionals (LC-wPBE, wB97, wB97X, wB97X-D, 
CAM-B3LYP) are found to outperform global hybrid 
functionals (B3LYP, PBE1PBE), pure density functional 
(LDA) and the exact-exchange (HF) method for excita- 
tion energy, oscillator strength and transition origin pre- 
diction, when compared with the EOM-CCSD results. It 
is also notable that the problematic description of Type- 
C, F and G excitations by B3LYP and PBE1PBE all 



concerns with N2-orbital involvement, which is HOMO- 
3, the deepest used standard-orbital from B3LYP KS- 
orbitals of Thy, and Type-C, G also concerns with S2- 
orbital, the LUMO+1, signifying that as (often higher- 
lying) electronic excitations are concerned with orbitals 
more far away from HOMO or LUMO, the correction of 
self-interaction error may be more important even if the 
system is as small as Thy. 
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TABLE SII. The standard-orbitals projection coefficients (coefficient square for orbital-fraction from backbone) of NT01(2)- 
H(E) of several SO-hosted excitations calculated by (a) EOM-CCSD and TDDFT with LC hybrid functionals, and (b) several 
other theoretical methods shown in Table II and Table SI. A(nm) and f column record absorption wavelengths and oscillator 
strengths, respectively, a column records <te for EOM-CCSD and <jm for other theoretical methods. In addition, there is data 
of the average a shown in the method column. 



(a) 



EOM- 
CCSD 

dT-0.029) 



TD- 
CAM- 
B3LYP 
(°"A/ : 
Thy-0.027; 
dT-0.048) 



TD- 
uB97X-D 

Thy-0.032; 
dT-0.088) 



TD- 
wB97X 

Thy-0.018; 
dT-0.035) 



TD- 
uB97 

Thy-0.028; 
dT-0.040) 



TD- 
LC-ujPBE 

Thy-0.029; 
dT-0.040) 



Thy 
N A(nm) 



f 



1 241.40 

2 218.80 

3 186.94 

4 180.36 

5 171.51 

6 160.85 
8 156.00 



0.0001 
0.2220 
0.0000 
0.0802 
0.2382 
0.0001 
0.0003 



1 249.48 

2 233.87 

3 195.27 

4 187.14 

5 180.42 

6 176.17 

7 166.72 
9 155.18 



0.0002 
0.1791 
0.0000 
0.0683 
0.1592 
0.0000 
0.0002 
0.3904 



1 249.72 0.0002 

2 234.05 0.1764 

3 194.84 0.0000 

4 188.03 0.0673 

5 180.66 0.1580 

6 177.27 0.0001 

7 167.82 0.0002 
9 155.06 0.4063 



1 242.13 

2 228.82 

3 190.78 

4 178.49 

5 174.43 

6 164.14 

7 158.49 
9 151.73 



0.0002 
0.2013 
0.0000 
0.0659 
0.1980 
0.0002 
0.0003 
0.3700 



1 238.63 

2 226.03 

3 189.05 

4 173.59 

5 170.20 

6 158.03 

7 153.85 
9 149.65 



0.0002 
0.2140 
0.0000 
0.0642 
0.2295 
0.0004 
0.0004 
0.3419 



1 241.89 

2 226.50 

3 191.81 

4 174.37 

5 171.41 

6 159.28 

7 155.34 

8 150.72 



0.0002 
0.2118 
0.0000 
0.0663 
0.2327 
0.0004 
0.0004 
0.3484 



NT01(2)-H 
PI Nl P2 N2 



-0.01 0.96 

0.99 0.01 

-0.01 0.59 

-0.01 0.01 

1.00 0.01 

0.00 0.27 

-0.01 0.83 



-0.01 
0.01 
-0.01 
1.00 
0.01 
-0.01 
0.00 



-0.23 
0.00 
0.80 
0.00 
0.00 
0.95 
-0.53 



-0.04 0.97 

0.99 0.04 

0.00 0.60 

0.08 0.00 

1.00 0.01 

0.00 0.28 

-0.01 0.78 

-0.08 0.00 



0.00 
-0.08 
0.00 
1.00 
0.08 
0.00 
0.00 
0.99 



-0.25 
0.00 
0.79 
0.00 
0.00 
0.95 
-0.62 
0.00 



-0.04 0.97 0.00 -0.25 

0.99 0.04 -0.08 0.00 

-0.02 0.61 0.00 0.78 

0.08 0.00 0.99 0.00 

0.99 0.01 0.09 0.00 

0.00 0.29 0.00 0.94 

0.00 0.76 0.00 -0.65 

-0.08 0.00 0.99 0.00 



-0.05 0.96 

1.00 0.05 

0.00 0.57 

0.05 0.00 

1.00 0.01 

0.00 0.26 

-0.01 0.83 

-0.05 0.00 



0.00 
-0.04 
0.00 
0.99 
0.03 
0.00 
0.00 
0.99 



-0.28 
0.00 
0.81 
0.00 
0.00 
0.94 
-0.55 
0.00 



-0.05 0.95 

1.00 0.05 

0.00 0.55 

0.08 0.00 

1.00 0.01 

0.00 0.24 

-0.01 0.85 

-0.03 0.00 



0.00 
-0.02 
0.00 
0.99 
-0.04 
0.00 
0.00 
0.99 



-0.29 
0.00 
0.83 
0.00 
0.00 
0.91 
-0.49 
0.00 



-0.04 0.95 

1.00 0.04 

0.00 0.54 

0.09 0.00 

1.00 0.01 

0.00 0.24 

-0.01 0.85 

-0.03 0.00 



0.00 
-0.02 
0.00 
0.99 
-0.04 
0.00 
0.00 
0.99 



-0.29 
0.00 
0.83 
0.00 
0.00 
0.91 
-0.49 
0.00 



NT01(2)-E 
SI S2 



0.94 0.28 

0.99 -0.01 

-0.36 0.91 

0.99 -0.01 

0.01 0.99 

0.94 0.31 

-0.27 0.92 



0.94 0.27 

1.00 -0.03 

-0.46 0.86 

1.00 -0.04 

0.08 1.00 

0.91 0.42 

-0.28 0.94 

-0.01 1.00 



0.95 0.26 

1.00 -0.03 

-0.49 0.85 

1.00 -0.05 

0.08 1.00 

0.90 0.43 

-0.25 0.95 

-0.02 1.00 



0.92 
1.00 
-0.35 
1.00 
0.06 
0.95 
-0.34 
-0.08 



0.31 
-0.05 
0.90 
0.06 
1.00 
0.30 
0.90 
1.00 



0.91 
1.00 
-0.31 
0.99 
0.06 
0.97 
-0.36 
-0.14 



0.33 
-0.06 
0.91 
0.12 
1.00 
0.22 
0.89 
0.99 



0.91 
1.00 
-0.31 
0.99 
0.06 
0.97 
-0.37 
-0.14 



0.34 
-0.06 
0.91 
0.13 
1.00 
0.22 
0.88 
0.99 



0.012 
0.028 
0.033 
0.028 
0.029 
0.033 
0.030 



0.013 
0.028 
0.042 
0.029 
0.031 
0.037 
0.042 



0.021 
0.022 
0.009 
0.027 
0.016 
0.006 
0.022 



0.027 
0.021 
0.021 
0.046 
0.021 
0.031 
0.030 



0.028 
0.019 
0.023 
0.050 
0.021 
0.031 
0.034 



dT 

N A(nm) 



f 



1 239.47 

2 225.95 

3 186.75 

4 179.74 

5 177.08 

6 164.55 

7 156.79 



0.0000 
0.3051 
0.0004 
0.1146 
0.2204 
0.0009 
0.0003 



1 247.05 

2 239.24 

3 194.94 

4 186.50 



0.0000 
0.2568 
0.0004 
0.1413 



5 184.44 0.1035 

6 180.74 0.0029 

7 169.05 0.0067 

8 166.90 0.0014 



1 247.18 0.0000 

2 239.44 0.2531 

3 194.60 0.0004 

4 187.12 0.1232 

5 184.84 0.1186 



6 181.44 

7 169.70 

8 167.85 



0.0024 
0.0064 
0.0018 



1 239.80 

2 233.96 

3 190.28 

4 179.36 

5 176.86 

6 168.19 

7 159.11 

8 156.42 



0.0000 
0.2820 
0.0003 
0.1967 
0.0825 
0.0015 
0.0010 
0.0086 



1 236.44 

2 231.07 

3 188.48 

4 175.51 

5 171.75 

6 162.13 

7 154.77 

8 150.79 



0.0000 
0.2966 
0.0002 
0.1967 
0.1108 
0.0016 
0.0013 
0.0208 



1 239.42 

2 231.30 

3 191.01 

4 176.69 

5 172.37 

6 163.08 

7 156.01 

8 151.61 



0.0000 
0.2938 
0.0002 
0.2022 
0.1125 
0.0014 
0.0010 
0.0317 



PI 



NT01(2)-H 
Nl P2 N2 



(B)' 



-0.02 
0.98 
-0.01 
-0.02 
0.98 
0.03 
-0.02 



0.96 
0.02 
0.69 
0.04 
0.02 
0.21 
0.84 



-0.03 -0.22 

0.02 0.00 

-0.05 0.66 

1.00 0.02 

0.02 0.00 

-0.05 0.79 

-0.02 -0.48 



0.01 
0.02 
0.05 
0.00 
0.02 
0.28 
0.04 



-0.01 
0.98 
0.04 
0.97 
-0.10 
0.16 
0.09 
-0.10 
0.00 



0.96 
0.01 
0.60 
-0.02 
-0.03 
0.00 
0.19 
0.09 
0.77 



0.00 -0.25 

-0.06 0.00 

0.00 0.76 

0.10 -0.03 

0.99 -0.01 

0.97 0.01 

0.00 0.73 

-0.03 0.61 

-0.01 -0.61 



0.01 
0.03 
0.05 
0.03 
(1.(11 
0.01 
0.36 
0.53 
0.02 



-0.01 
0.98 
0.04 
0.89 
-0.41 
0.46 
0.86 
0.09 
-0.10 
0.00 



0.97 
0.01 
0.60 
-0.03 
-0.01 
-0.01 
-0.02 
0.21 
0.05 
0.75 



0.00 -0.25 

-0.06 0.01 

0.00 0.75 

0.41 -0.03 

0.91 0.01 

0.86 0.00 

-0.48 -0.02 

0.01 0.75 

-0.04 0.58 

-0.01 -0.64 



0.00 
0.03 
0.05 
0.03 
0.01 
0.03 
0.02 
0.31 
0.56 
0.02 



-0.01 
0.98 
0.03 
0.98 
0.02 
0.05 
0.00 
-0.08 



0.96 
0.01 
0.57 
-0.01 
-0.01 
0.20 
0.78 
0.28 



0.00 -0.28 

-0.02 0.00 

-0.01 0.79 

0.04 -0.02 

0.99 0.00 

0.00 0.74 

0.02 -0.55 

-0.05 0.50 



0.00 
0.03 
0.03 
0.03 
0.00 
0.29 
0.03 
0.54 



-0.01 
0.98 
0.03 
0.97 
-0.03 
0.04 
-0.01 
-0.05 



0.95 
0.02 
0.55 
-0.01 
-0.01 
0.20 
0.79 
0.30 



0.00 -0.29 

-0.01 0.00 

-0.01 0.81 

0.05 -0.02 

0.99 0.00 

0.01 0.72 

0.02 -0.50 

-0.09 0.49 



0.00 
0.03 
0.03 
0.03 
0.00 
0.29 
0.03 
0.51 



-0.01 
0.98 
0.03 
0.97 
-0.02 
0.04 
-0.01 
-0.04 



0.95 
0.02 
0.54 
0.01 
-0.01 
0.20 
0.81 
0.29 



0.00 -0.29 

-0.01 0.00 

-0.01 0.82 

0.05 -0.02 

0.99 0.00 

0.01 0.72 

0.02 -0.49 

-0.12 0.49 



0.00 
0.03 
0.03 
0.03 
0.00 
0.28 
0.03 
0.51 



NT01(2)-E 
SI S2 



0.95 0.20 

0.98 0.00 

-0.33 0.92 

0.98 0.00 

-0.01 0.97 

0.96 0.16 

-0.35 0.90 



0.94 0.28 

0.99 -0.02 

-0.47 0.85 

-0.09 0.98 

0.98 0.08 

0.97 0.16 

0.91 0.35 

0.98 -0.01 

-0.26 0.94 



0.94 0.27 

0.99 -0.02 

-0.50 0.84 

-0.40 0.90 

0.90 0.39 

0.86 0.47 

-0.49 0.86 

0.90 0.38 

0.98 -0.06 

-0.23 0.95 



0.92 0.32 

0.99 -0.03 

-0.34 0.90 

0.04 0.98 

0.99 0.08 

0.95 0.24 

-0.48 0.84 

0.97 0.06 



0.91 0.34 

0.99 -0.04 

-0.31 0.91 

0.07 0.98 

0.98 0.14 

0.96 0.19 

-0.51 0.81 

0.97 0.04 



0.90 0.34 

0.99 -0.04 

-0.30 0.91 

0.06 0.98 

0.98 0.15 

0.96 0.19 

-0.50 0.82 

0.97 0.03 



0.024 
0.006 
0.052 
0.012 
0.011 
0.067 
0.029 



0.026 
0.024 
0.063 
0.036 

0.071 
0.064 

0.052 



0.024 
0.024 
0.069 
0.163 



0.065 



0.041 
0.015 
0.054 
0.019 
0.030 
0.032 
0.051 



0.048 
0.015 
0.062 
0.027 
0.043 
0.028 
0.056 



0.048 
0.015 
0.066 
0.023 
0.046 
0.028 
0.051 



SG 



Thy 
N A(nm) 



f 



PI 



NT01(2)-H 
Nl P2 



X2 



NT01(2)-E 
SI S2 



dT 

N X(nm) f 



PI 



NT01(2)-H 
Nl P2 N2 (B) s 



NT01(2)-E 
SI S2 



CIS 

{"M ■ 
Thy-0.138; 
dT-0.080) 



197.73 0.0077 
195.68 0.4332 
161.65 0.0000 



4 147.53 0.1508 

6 141.31 0.4179 

8 126.65 0.2298 

12 116.52 0.0004 

15 110.63 0.0001 



-0.13 
0.98 
-0.01 
0.80 
-0.59 
-0.61 
0.13 
0.00 
-0.01 
0.00 
-0.02 



0.93 -0.01 

0.12 0.02 

0.47 -0.01 

0.02 0.58 

0.00 0.79 

0.00 0.78 

0.01 0.97 

0.59 -0.01 

0.63 0.00 

0.03 0.00 

0.69 0.00 



-0.33 
-0.04 
0.87 
0.00 
0.00 
0.01 
0.00 
0.61 
-0.66 
0.63 
-0.40 



0.86 0.40 

0.99 -0.07 

-0.23 0.91 

0.56 0.80 

0.79 -0.58 

0.62 0.77 

-0.38 0.92 

0.98 0.01 

0.05 0.93 

0.98 0.17 

-0.19 0.92 



0.062 
0.038 
0.055 
0.242 



1 198.84 0.5502 

2 195.34 0.0001 

3 160.17 0.0005 

4 149.55 0.2503 



141.28 0.2894 
126.31 0.2175 



15 115.65 0.0033 



0.98 0.01 0.03 0.00 0.02 

-0.02 0.94 -0.01 -0.32 0.00 

0.03 0.48 -0.02 0.86 0.02 

0.98 0.00 0.07 -0.03 0.03 

-0.21 0.00 0.96 0.01 0.00 

0.10 0.03 0.97 -0.01 0.00 

0.04 -0.46 0.01 -0.56 0.07 

0.90 0.09 0.14 0.02 0.13 



0.98 -0.06 

0.85 0.41 

0.23 -0.91 

-0.10 -0.98 

0.85 0.49 

0.41 -0.90 

0.96 0.03 

-0.08 0.05 



0.018 
0.074 
0.090 
0.037 
0.157 



TD-HF 

(o"M : 
Thy-0.137; 
dT-0.106) 



1 206.52 0.3427 

2 202.72 0.0021 

3 164.91 0.0000 

4 152.36 0.1127 

5 144.72 0.3622 
8 129.32 0.1426 
15 111.24 0.0004 
22 100.71 0.0099 



0.98 
0.05 
-0.01 
0.65 
-0.75 
0.66 
0.18 
0.00 
0.61 
-0.02 



-0.04 0.04 

0.94 0.00 

0.47 -0.01 

0.01 0.73 

-0.01 0.64 

0.00 -0.74 

0.03 0.95 

0.08 0.00 

0.06 0.05 

0.88 -0.01 



0.02 
-0.33 
0.87 
0.00 
0.00 
0.00 
-0.01 
0.65 
0.00 
0.44 



0.99 -0.07 

0.85 0.41 

-0.23 0.91 

0.70 0.67 

0.67 -0.73 

0.53 0.84 

-0.41 0.90 

0.97 0.21 

-0.13 0.70 

0.01 -0.04 



0.025 
0.057 
0.055 
0.297 



1 209.49 0.4576 

2 200.20 0.0000 

3 163.31 0.0005 

4 153.99 0.1954 

5 144.95 0.2422 
7 128.95 0.1427 
15 116.43 0.0077 

21 110.25 0.0317 



0.97 0.01 0.04 0.00 0.01 

-0.02 0.94 -0.01 -0.32 0.00 

0.04 0.48 -0.02 0.86 0.02 

0.97 0.00 0.09 -0.04 0.02 

-0.31 0.01 0.93 0.02 0.00 

0.14 0.04 0.96 -0.01 0.00 

0.50 -0.30 0.07 -0.39 0.28 

0.70 0.32 0.08 0.34 0.04 

0.96 0.03 0.01 0.00 0.02 



0.98 -0.06 

0.84 0.41 

0.23 -0.91 

-0.11 -0.97 

0.79 0.59 

0.43 -0.89 

0.69 0.06 

0.64 -0.02 

0.87 0.25 



0.019 
0.075 
0.090 
0.042 
0.199 



TD- 
B3LYP 

( CT M : 
Thy-0.116; 
dT-0.110) 



1 270.56 0.0001 

2 246.47 0.1268 

3 216.14 0.0000 

4 208.21 0.0574 

5 200.48 0.0001 

6 192.17 0.1223 

7 184.80 0.0002 
9 162.81 0.3352 



-0.04 
0.97 
0.00 
0.23 
-0.05 
0.99 
0.03 
-0.16 



0.99 0.00 

0.04 -0.21 

0.25 0.00 

0.00 0.97 

0.91 0.00 

0.04 0.15 

0.38 0.00 

0.00 0.98 



-0.15 
0.00 
0.96 
0.00 
0.40 
0.05 
-0.92 
0.00 



0.99 0.15 

1.00 0.00 

0.97 -0.24 

0.97 -0.22 

0.10 0.99 

0.14 0.99 

-0.03 1.00 

0.01 1.00 



0.048 
0.064 
0.159 
0.093 
0.200 
0.058 
0.187 



1 


268.35 


0.0001 


2 


253.13 


0.1897 


3 


220.32 


0.0042 


4 


207.29 


0.0591 


5 


205.43 


0.0003 



6 200.39 0.0014 



7 197.43 0.1530 
9 186.15 0.0003 



-0.03 
0.94 
0.15 
0.26 
-0.16 
0.01 
0.08 
0.21 
0.96 
-0.04 



0.98 
0.02 
0.11 
-0.03 
0.53 
0.69 
0.67 
0.62 
-0.04 
0.42 



0.00 
-0.16 
-0.02 
0.96 
0.12 
-0.15 
0.11 
-0.07 
0.18 
0.00 



-0.15 0.01 

-0.03 0.07 

0.74 0.37 

-0.04 0.01 

0.53 0.39 

0.12 0.45 

0.68 0.07 

-0.43 0.35 

-0.04 0.04 

-0.87 0.06 



0.97 
0.99 
0.98 
0.95 
-0.87 
0.46 
0.60 
-0.78 
0.19 
-0.07 



0.16 
0.01 
-0.10 
-0.28 
0.47 
0.86 
0.77 
0.61 
0.95 
0.98 



0.026 
0.054 
0.089 
0.118 
0.221 



0.077 
0.186 



TD- 
PBE1PBE 

(ctm : 
Thy-0.113; 
dT-0.096) 



1 263.26 0.0001 

2 241.45 0.1414 

3 207.93 0.0000 

4 201.76 0.0589 

5 193.99 0.0002 

6 187.83 0.1333 

7 179.67 0.0002 
9 159.34 0.4077 



-0.04 0.98 0.00 -0.18 

0.98 0.04 -0.17 0.00 

-0.01 0.37 0.00 0.92 

0.19 0.00 0.98 0.00 

-0.05 0.79 0.00 0.61 

0.99 0.04 0.13 0.05 

0.02 0.49 0.00 -0.87 

-0.13 0.00 0.99 0.00 



0.98 0.18 

1.00 -0.01 

0.91 -0.41 

0.98 -0.18 

0.32 0.95 

0.12 0.99 

-0.08 1.00 

0.02 1.00 



0.036 
0.053 
0.210 
0.076 
0.212 
0.050 
0.151 



1 261.13 0.0001 

2 247.82 0.2086 

3 211.64 0.0033 

4 200.76 0.0643 

5 199.04 0.0001 

6 193.41 0.0911 

7 192.11 0.0688 
9 180.60 0.0003 



-0.03 
0.95 
0.13 
0.22 
-0.11 
-0.10 
0.79 
-0.40 
-0.41 
0.63 
-0.02 



0.98 
0.02 
0.17 
-0.03 
0.82 
-0.13 
0.38 
-0.17 
0.29 
-0.47 
0.53 



0.00 
-0.13 
-0.01 
0.97 
0.02 
0.19 
0.23 
0.61 
-0.28 
-0.15 
0.01 



-0.18 0.01 

-0.02 0.06 

0.74 0.37 

-0.03 0.01 

0.54 0.03 

0.02 0.89 

0.35 0.05 

0.52 0.14 

0.69 0.20 

0.17 0.31 

-0.82 0.04 



0.97 
0.99 
0.97 
0.96 
-0.37 
0.91 
0.32 
0.93 
0.87 
-0.46 
-0.11 



0.19 
0.00 
-0.17 
-0.24 
0.92 
0.35 
0.92 
-0.33 
0.45 
0.87 
0.98 



0.017 
0.045 
0.102 
0.102 
0.063 



TD- 
LDA 

Thy-0.201) 



1 316.41 0.0000 

2 270.51 0.0662 

3 267.18 0.0003 

4 236.19 0.0011 

5 235.14 0.0507 

6 208.80 0.0944 

7 203.29 0.0001 
11 173.87 0.2857 



-0.03 1.00 0.00 0.05 

0.90 0.03 -0.40 0.05 

-0.05 -0.05 0.01 0.99 

-0.08 0.99 -0.03 -0.07 

0.43 0.06 0.90 0.01 

0.97 0.04 0.20 0.02 

-0.02 0.06 -0.01 0.99 

-0.19 0.00 0.97 0.01 

0.34 -0.01 0.16 0.00 



1.00 -0.04 

0.99 0.08 

1.00 -0.01 

0.04 1.00 

0.91 -0.40 

0.22 0.97 

0.00 1.00 

-0.05 1.00 

1.00 0.05 



0.125 
0.125 
0.133 
0.170 
0.174 
0.083 
0.195 



1 316.27 0.0006 

2 291.65 0.0413 

3 271.79 0.0622 

5 254.91 0.0125 

6 239.34 0.0006 

7 235.80 0.0456 

9 220.59 0.0679 

10 214.82 0.0548 



-0.07 
0.50 
0.67 
-0.33 
-0.08 
0.40 
0.42 
0.78 



0.95 
0.27 
-0.08 
0.07 
0.97 
0.05 
0.15 
-0.05 



0.02 
-0.16 
-0.25 
0.14 
-0.02 
0.86 
0.03 
0.15 



0.07 0.07 

-0.28 0.53 

0.59 0.11 

0.73 0.34 

-0.05 0.04 

-0.03 0.09 

-0.33 0.65 

0.35 0.22 



0.99 
0.99 
0.98 
0.99 
0.01 
0.90 
0.02 
0.11 



-0.02 
0.01 
0.05 
0.03 
0.98 
-0.41 
0.97 
0.93 



